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Multipolar interactions are derived in the orbitally degenerate Anderson lattice with a spher- 
ical Fermi surface and one conduction electron per unit cell of the simple cubic lattice. As the 
crystalline-electric-field (CEF) ground state of / configuration, the four-fold degenerate Y% is 
mainly studied. Intersite interactions up to a sufficiently distant pair are Fourier transformed 
to the wavenumber space. For the Fg case, quadrupolar and octupolar interactions favor the 
staggered order with q — (1/2, 1/2, 1/2), while the dipolar interaction favors an incommensu- 
rate magnetic structure with a long modulation period. The latter is due to a Kohn anomaly 
which is found to be sharply peaked for interaction channels with large angular momenta. 
Implications of results are discussed for multipole orders in CeBs, and the incommensurate 
magnetic structure in a quasi-cubic system CeB2C2- 
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1. Introduction 

It is a long-standing problem in magnetism to identify 
the microscopic origin of complicated magnetic struc- 
tures in solids. In some cases, a complicated structure 
results from coupling between dipoles and orbital de- 
grees of freedom, or the quadrupolc moment. Accord- 
ingly ordering of multipole moments higher than dipoles 
has attracted much attention. In CeBg for example, 
a peculiar magnetic structure is realized on the back- 
ground of a quadrupole (or orbital) order. 1 ) The resul- 
tant wave number q = (1/4, ±1/4, 1/2) of the magnetic 
moment is closely related to the quadrupole order with 
q = (1/2, 1/2, 1/2). Thus it is tempting to interpret the 
magnetic pattern in terms of the background quadrupole 
patterns together with the strong spin-orbit interaction. 

However, a magnetic order with q = (1/4, ±1/4, 1/2) 
has also been found in other RBg systems including R 
= Gd. 2 - 1 Since the trivalent Gd with 4/ 7 configuration 
does not have the orbital degrees of freedom, the char- 
acteristic q should originate from structure of conduc- 
tion bands. A model has been proposed to explain the 
stability of q = (1/4, ±1/4, 1/2). 3 ) The model empha- 
sizes the interpocket polarization among the three nearly 
spherical Fermi surfaces in the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction. 

On the other hand, related compounds RB2C2 (R = 
rare earths), which crystallize in the tetragonal LaB2C2- 
type structure(P4/m&?7i), 4 ) often exhibit an incommen- 
surate magnetic structure. 5 - ) For example, CeB2C2 un- 
dergoes a magnetic ordering at Xn = 7.3K, and a suc- 
cessive transition at T t = 6.5K. 6 ) Below T t , a long- 
period magnetic structure with its propagation vector 
k = (6, 5, S'), where 5 = 0.16 and 6' = 0.10, is found from 
neutron diffraction. 5 ) Thus propagation vectors seem to 
be characterized by conduction bands near the Fermi 
level. 

RB2C2 can be considered as a quasi-cubic system be- 
cause the distance a/y/2 between R ions in the c-plane is 
almost the same as the distance c along the c-axis. 4 - 1 The 



underlying atomic composition is closely related to RBg 
systems. Namely, if one removes two boron atoms from 
each B6 cluster, and further replaces two boron atoms 
by two carbon atoms to make a B2C2 cluster, one ob- 
tains the RB2C2 structure. According to band-structure 
calculation for LaB2C2, 7 - ) a connected ellipsoidal Fermi 
surface is formed around the Z point in the simple tetrag- 
onal Brillouin zone (BZ). This Fermi surface is partly 
confirmed by de Haas-van Alphen (dHvA) experiment. 8 ) 
Qualitatively speaking, the ellipsoidal shape is a result 
of growth from one of three equivalent Fermi surfaces in 
RBg. The other two pieces shrink and disappear. 

It should be clear from the examples above that the 
interplay of energy-band structure and the 4/ shell struc- 
ture deserves detailed inspection. The purpose of this pa- 
per is to derive intersite interactions of multipoles from 
microscopic calculation for a simple model which still 
keeps the realistic 4/ shell structure. The RKKY inter- 
action should largely be determined by hybridization in 
the case of CeB 6 and related compounds. In a previ- 
ous paper, we calculated the RKKY interaction for the 
nearest and the next nearest neighbors taking the An- 
derson lattice with orbital degeneracy. The conduction- 
band states are taken as plane waves and the hybridiza- 
tion as spherically symmetric. 9 ) In this paper we extend 
the calculation to many different pairs to obtain the in- 
tersite interaction in the wavenumber space. We again 
take the spherical Fermi surface in the simple cubic BZ 
(SCBZ) with the lattice constant a q = a/y/2. The result 
can be mapped onto the the quasi-cubic BZ (QCBZ) of 
the bace-centred tetragonal lattice with lattice constants 
a and c — a q . We consider mainly the r§ quartet of total 
angular momentum j — 5/2 as the lowest CEF states. 
This quartet includes the m = ±1/2 doublet where m is 
the azimuthal quantum number of the total angular mo- 
mentum along the c axis. The doublet is likely to be the 
CEF ground state in the tetragonal CeB2C2. 10 - ) Hence we 
shall also derive corresponding results for the m — ±1/2 
doublet, and for another doublet IV 

This paper is organized as follows. In §2 we briefly 
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review the effective intersite Hamiltonian derived previ- 
ously. 9 ^ In §3 multipolar couplings are calculated numer- 
ically for the simple cubic lattice, and are transformed 
to wavenumber space. Relation between possible order- 
ing structures and the Fermi surface is discussed with 
emphasis on the Kohn anomaly. We discuss implications 
of results on actual materials in §4. Appendix gives some 
analytic results about the polarization functions with fi- 
nite angular momentum. 

2. Intersite Interactions from Real to Wavenum- 
ber Spaces 

We start from the spherically symmetric form of hy- 
bridization: 



an excited state £. If we take £ as the 4/° state or a 4/ 2 
state with orbital angular momentum L, spin S and total 
angular momentum J, A(£; j t , j M ,P) is given by 

A(/ ;i t ,i„,P)=-Oi||/jj|0)(j 1 ||/jj|0)* 



_lV'i+i«- p 



x (-!> 



h h p 
it j u o 



(3) 

A(/ 2 (L5)J;j t ,i„,P)=((L5)J||4||i 1 )*((L5)J||/tj|j 1 ) 

x (-i)h+j t +j+2P fji h P\ 
lit j u J) 

(4) 

From eq.(2), we obtain the intersite interaction pro- 
jected by V to low-lying CEF states: 



ihnur 



(1) ^VHi/P 



where N is the number of lattice sites, V is the strength 
of hybridization, and Y; m (Q) is the spherical harmonics 
with solid angle O. A 4/ electron is created by fm a {i) at 
site R{ with z-component m of the orbital angular mo- 
mentum I = 3 and spin a. A conduction electron with 
wavenumber k and spin a is annihilated by c^ a . Each 
Ce 3+ ion has one 4/ electron, and the spin-orbit cou- 
pling leads to the total angular momentum j\ = 5/2 in 
the ground state. By using the fourth-order perturbation 
theory with respect to i?hyb, we obtain the multipolar in- 
tersite interaction between H and H 7 -: 9 ) 
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where X a (i) is a multipole operator at site H, 9 ) 
D a ^(^,^';Rij) is the dimensionless multipolar coupling 
strength. In the case of r 8 state, a and f} run from 1 to 
15. Their Fourier transforms are defined as follows: 
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where fep is the Fermi wavenumber, and Ri 
We have introduced E a — h 2 /2moal with mo being the 
effective mass of conduction electrons. In eq.(2) Qi^ix) 
is the range function tabulated in ref.9, -E e xc(£) is the 
excitation energy to a state £. The shorthand notation 
■ • • ,3n] = (2ji + 1) • • • (2j n + 1) is used, and round 
or curly brackets denote Wigner's 3n-j symbols. We have 
introduced T p (i) as the set of rank-P irreducible-tensor 
operators at site H with (ji||T p ||ji) its reduced matrix 
element, and {■ ■ ■ }§ denotes a tensor component of rank 
zero, i.e., a scalar. n ) The factor A(£;j t ,j„,P) includes 
information of wave functions of the ground state and 



N ■ 
q 

(7) 

With the time-reversal symmetry, multipole operators 
are classified into electric (even) and magnetic (odd) mul- 
tipole operators. In the intersite interaction, the opera- 
tors with different time-reversal symmetry do not cou- 
ple. Therefore we obtain a block-diagonalized matrix for 
D ai3 (£,,£,';q) where submatrices £>„(£,£';<?) with either 
even (v = g) or odd {y = u) symmetry under the time 
reversal can be defined. Then we diagonalize the subma- 
trices for each q. We note that comparison of relative 
strength of resultant multipolar couplings is physically 
relevant because all A Q (i)'s are normalized in the same 
way. In the next section, we diagonalize £>„(£,£'; <z) nu- 
merically and discuss possible ordering structures. 

3. Numerical Results of Intersite Interactions 

To obtain the wavenumber dependence of intersite in- 
teractions, we first calculate D v (^, Roj) numerically 
for a sufficiently large number of sites Rj . As for the in- 
termediate states, both 4/° and 4/ 2 Hund's-rule ground 
states are considered, and they are referred to f° and 
/h respectively. We define three cases depending on the 
intermediate states at Ho and H,: (a) f° intermediate 
states on both sites, (b) intermediate states on both 
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sites, and 

(c) f° on one site and on the other sites. We ex- 
clude D u (£, R = 0) in Fourier transform because this 
belongs to the on-site interaction. After taking Fourier 
transform, we diagonalize D v (^,^';q) for each q. We 
have compared numerical results by varying the number 
N of neighbors as N = 32 3 , 64 3 and 128 3 , and found 
a reasonable convergence. The remaining difference is 
consistent with the asymptotic range dependence of the 
RKKY interaction given by R~ 3 . We have also checked 
numerically that D v {£,, C: Q) is Hermitian within relative 
error of order 10~ 5 . In the following analysis we use the 
results with N = 64 3 . 

3. 1 Quadrupolar Interaction 

We begin with the quadrupolar interactions that 
are even under time reversal in the r§ manifold. Fig- 
ure 1 shows five eigenvalues of quadrupolar couplings 
D g (t,£;q) in the SCBZ with k F a q = (37T 2 ) 1 / 3 , w hi c h 
corresponds to one conduction electron per cubic unit 
cell. The absolute minimum gives the most stable prop- 
agation vector and symmetry of the electronic order. At 
high symmetry points such as T and R, eigenvalues are 
characterized by irreducible representations of multipoles 
such as T^g and T§ g . This is consistent with the symme- 
try analysis of ref.13. For all cases (a), (b) and (c), the 
minimum occurs at the R point and the most favorable 
symmetry is T^g. In the case (a), the eigenvalues of T^ g 
are comparable at the R point. On the other hand, T 3g 
dominates the stability in the cases (b) and (c). In the 
previous analysis with only the nearest-neighbor sites, 
the r 5s order was found more stable than the r 3g for the 
case (a). 9 ) The difference from the present result means 
that interactions between third nearest neighbors or be- 
tween further distant sites pile up to change the most 
stable order. It remains to identify the reason why the 
actual CeBg prefers the T§ g order. Certainly the actual 
Fermi surface with three identical pieces should influence 
the selection. 

We note the sharp extrema near the T point in Fig. 
1, which arise due to the Kohn anomaly. As we have 
2fcp / '(2ir I 'a q ) ~ 0.98, the location of the Kohn anomaly 
is near the T point in the reduced zone scheme. It is well 
known that the Kohn function (or the static limit of the 
Lindhard function) does not show a peak at 2fcp, and 
then the Kohn anomaly appears only as its divergent 
derivative. With orbital angular momenta as in eq.(l), 
however, some generalized polarization functions show a 
hump or a dip around 2fcp as discussed later. 

3.2 Magnetic Interaction 

Figure 2 shows eigenvalues of magnetic multipolar cou- 
plings £>„(£, in the SCBZ with k F a q = (3tt 2 ) 1/3 - 
The minimum gives the propagation vector and symme- 
try of the favorable magnetic order. Ten magnetic multi- 
pole operators in the r§ quartet are characterized by ir- 
reducible representations T 2u , r 4 „i, r 4 „2 and T 5u . While 
the multipolar operators of T^i or r4„2 symmetry both 
have magnetic dipole moment, T 2 u and T^ u do not con- 
tain the dipole character. The eigenvalues of r 4u2 and 
r 5u are degenerate at T and R points for the case (a). 




m x r r x r MR 
Wavenumber q 




f* M M ?. r xr mr 
£ Wavenumber q 



Fig. 1. Wavenumber dependence of five eigenvalues of quadrupo- 
lar interaction £? 9 (£,£';q). Results are shown for cases (a) f° 
intermediate states on both sites, (b) intermediate states on 
both sites, and (c) f° on one site and /jt on the other. 



For the cases (b) and (c), the magnetic interactions are 
mixtures corresponding to T^i and symmetries. 

For the case (a), the minimum appears near the T 
point, and its symmetry is close to T4 ul . The wavenum- 
ber is q K = (8 ,6 ,5 ) with S = 0.03 in units of 2n/a q , 
which is close to 1 — 2fcp / '(2n / 'a g ) ~ 0.02. Comparing the 
minima in Figs. 1 and 2, we find that the corresponding 
magnetic order is more stable than the quadrupole order. 
This means that the incommensurate magnetic order is 
favored with the present set of parameters. We emphasize 
that the Kohn anomaly is responsible for the stability. 

With orbital angular momenta as in eq.(l), generalized 
polarization functions are given by 

£ 5 k - k > J {tk) ~ f{ek,) Y hmi (k)Y hm2 (k ) 
— ' efc — ey 



4 J. Phys. Soc. Jpn. 



Full Paper 



Gen'ya Sakurai and Yoshio KuRAMOTO 



_L ( k F a q \ ^ r h+h+i : 
£ F V J 



T4u2, T5u 



x / dfiii* 8 (fi)^ iroi (fi)r Iam2 (n), 



(8) 



where /(e) is the Fermi distribution function, is the 
energy dispersion of free electrons, and li,l 2 are even 
integers with < I1J2 < 6. Note that terms with l 3 7^ 
give anisotropic functions in the wavenumber space. 
In addition, angular momentum in eq.(8) makes the In- 
dependence of the polarization function rather different 
from the ordinary Kohn function. For understanding the 
characteristic behaviour for finite l\ and fo, we consider 
the hypothetical RKKY interaction angular momentum 
/ ^> 1 and with m = 0. The corresponding polarization 
function is given by 

/(<*)"/(**'), 




-0.02 
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where Yi (k) corresponds to Y 3m (k) in eq.(l). The sig- 
nificant contribution in summation over k and k comes 
only from directions around the z-axis. This is because 
the spherical harmonics with I ^> 1 behaves asymptoti- 
cally as 12 ) 
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(10) 

and is sharply peaked at = 0, 7r. Then the polariza- 
tion function acquires a one-dimensional character with 
I 3> 1 and is strongly peaked at q = (0,0, ±2fcp). This 
anisotropic behaviour should be reflected in the terms 
with I3 7^ in eq.(8) around \q\ — 2kp. Figure 3 
shows radial parts of generalized polarization functions 
X l ili 2 {p) with p = \q\/kF- Except for Xoo(p)' wri ich is the 
Kohn function, these functions have a hump or a dip 
at |q| ~ 2kp. They can accumulate to give a sharp ex- 
trcmum near \q\ = 2fcp in the intersite interaction. 

The wavenumber for the Kohn anomaly in the SCBZ 
is given by Sk = 1— 2k? / {2-k / a q ) in units of the reciprocal 
lattice. To be precise, Sk differs to some extent from 8 
obtained numerically, because each extremum in (p) 
is off from p — 2 in general. We note that the magnetic 
interaction takes another minimum at the M points (1/2, 
1/2, 0), (0, 1/2, 1/2) and (1/2, 0, 1/2), and its symmetry 
is close to T4„2 ■ The anisotropy of the exchange is small 
in the case (a), but is large in cases (b) and (c). This 
is consistent with the results obtained in the previous 
paper. 9 ) 

For cases (b) and (c), the minimum occurs at the R 
point, and its symmetry is T5 U . This means that the an- 
tifcrro octupolar ordering is stabilized. The wavenumber 
q K also gives a local minimum, but the anomaly is weaker 
than the case (a) . In the case (c) , the dipole interaction is 
insignificant for all wavenumbers. This is consistent with 
the smallness of dipolar couplings in the nearest and next 
nearest neighbors as derived in the previous paper. 9 ) 



< 7T %4 




m x r 



R x r 
Wavenumber q 



M R 



Fig. 2. Wavenumber dependence of ten eigenvalues of magnetic 
multipolar interaction D u (£, <f)- Results are shown for cases 
(a) f° intermediate states on both sites, (b) intermediate 
states on both sites, and (c) f° on one site and on the other. 




Fig. 3. Radial parts of generalized polarization functions X;J; 2 (p) 
with p = |q|/feF. 
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3.3 Other Choices of Lowest CEF States 

So far multipolar interactions are calculated with T% 
quartet as the lowest CEF states. In this subsection, we 
choose other states as the lowest CEF states: (i) j\ = 
5/2, m = ±1/2 doublet and (ii)r 7 doublet of ji = 5/2. 
Here only the 4/° is considered as an intermediate state. 

We first take the dimensionless wavenumber kpa q = 
0.987T, which corresponds to one conduction electron per 
cubic unit cell. Figure 4(a) shows magnetic dipolar in- 
teraction for the m = ±1/2 doublet. The wavefunction 
has uniaxial anisotropy along [001] direction. Hence the 
results are shown for a BZ accounting for the correspond- 
ing anisotropy. With 4/° intermediate states, the dipo- 
lar interaction between the Kramers doublets is a scalar. 
Around the T point, the Kohn anomaly appears and gives 
the minimum. The preferred wavenumber is the same as 
that for the Tg quartet as the lowest CEF states. We 
note that the Kohn anomaly is conspicuous along [001], 
while almost invisible along [100]. Such strong anisotropy 
comes from the highly elongated shape of the wave func- 
tion along [001]. 

Figure 4(b) shows magnetic dipolar interaction for the 
Tt doublet where the Kohn anomaly is not obvious. How- 
ever, the derivative of D u (f° , f°;q) with respect to q 

seems to diverge at (5^\ 4^, 4^) with 4^ - °- 43 - This 
is interpreted as the three-dimensional Kohn anomaly. 
Namely, taking into account the reciprocal lattice vector 
along [111], we obtain the Kohn anomaly at S' K = 1 — 
kpa q /\/3iT. Here the dimensionless wavenumber k-pa q is 
0.57-^71" and gives S' K ~ 0.43. Three more anomalies can 
be found at (5^,5^,0) with 5^ ~ 0.32, (4 3) , Sf\ 1/2) 
with 4 3) ~ 0.20, and (6^,6^,1/2) with 6 { 7 i] ~ 0.39. 
These wave numbers are interpreted as intercepts of 
spheres with its radius 2fcp at high symmetry axes in 
the irreducible BZ. The spheres are centered at recipro- 
cal lattice points [110], [100] and [111], respectively, 

In order to study the effect of band structures in the 
simplest way, we now change the Fermi wavenumber. 
Namely we take k?a q = 0.86n ~ 0.497\/37r, and show 
the results in Fig. 5. For comparison with the main case 
of kpa q = 0.987T, we also show in Fig. 5(a) the magnetic 
multipolar interaction with the T$ quartet. Figure 5(b) 
shows the magnetic dipolar interaction for the m = ±1/2 
doublet. The Kohn anomaly is found at (0,0, <5 ) with 
5 = 0.17, which is close to 1— kpaq/ir = 0.14. In the case 
of T 7 doublet shown in Fig. 5(c), on the other hand, Kohn 
anomalies in the magnetic dipolar interaction arc found 
around the R point. The minimum is at (5^,5^,5^) 

with 4 5) - 0.48. Since 2k F ~ 0.994\/37r/ a g is almost the 
same as the distance between T and R points, the Kohn 
anomaly can arise near the R point. From the large dif- 
ference between the results for each doublet, we see that 
the anisotropy of the lowest CEF states sensitively in- 
fluences the wavenumber dependence of the interaction, 
and the character of the Kohn anomaly. 

4. Discussion 

4-1 Comparison with real-space results 

We now compare between the present result and the 
previous one 9 ) which takes into account only the nearest 




-0.04 J 1 1 1 1 1 

MXT R XT MR 

Wavenumber q 



Fig. 4. Wavenumber dependence of magnetic interaction 
D u (f ,f°;q) with k F a q = 0.98n. Results arc shown for the 
lowest CEF states (a) j 1 = 5/2, m = ±1/2 doublet, and (b) T 7 
doublet. Arrows indicate the location of the Kohn anomaly. 



and the next-nearest neighbor interactions. The previous 
study in addition neglects possible mixing of different 
point-group representations in the intersite interaction. 
As clarified by Sakai et a/., 13 ) the mixing becomes min- 
imum at high symmetry points such as R in the SCBZ. 
Thus the result can be compared with previous study if 
we look at only T and R. However T^i and T4 U 2 repre- 
sentations still mix at T and R points as seen in Fig 2(b) 
and (c). Concerning the time-reversal odd interactions, 
occurrence of ferromagnetic dipolar interaction and an- 
tifcrromagnctic T 5u octupolar interaction is consistent 
with the previous result. In the present case, however, it 
is difficult to judge whether the interactions have a sig- 
nificant anisotropic, or "pseudo-dipolar" character. If we 
look at only the R point in Fig. 2(b) and (c), eigenvalues 
of appear as triply degenerate. The same degener- 
acy also applies to the eigenvalues of T^ u representations. 
The pseudo-dipole character can only be seen by gather- 
ing information at various wavenumbers. 

On the other hand, the implication for quadrupole or- 
der is different between the previous result 9 ^ and the 
present one. Namely we found in the previous study 
that the nearest-neighbor interaction prefers staggered 
quadrupole order with T 5g symmetry to T 3g , because 
the interaction between T 5g moments is antifcrro-type 
but that between T 3g moments includes both ferro- and 
antiferro-type couplings. In the present study, it is found 
that many distant pairs interact and make r 3g more fa- 
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Fig. 5. Wavenumber dependence of magnetic interaction 
f>„(f°,f°;q) with k F a q = 0.86vr. Results are shown for the 
lowest CEF states, (c) Tg quartet, (b) ji = 5/2, m = ±1/2 
doublet, and (c) T-j doublet. Arrows indicate the location of the 
Kohn anomaly. 



vorablc, although the R point remains the most favorable 
wavenumber. The relative stability of T^ g against T<=, g is 
fragile if f° is the dominant intermediate state. As the 
intermediate states become important, T^g becomes 
more stable. It is the effect of many distant neighbors 
that makes the difference from previous results. Thus 
there is no inconsistency between both calculations. 

4-2 Relevance to real Ce materials such as CeB 6 and 
CeB 2 C 2 

In CeB 6 , it is known that the quadrupole order sets 
in with the r 5g symmetry. It is not known, however, to 
what extent it is stable against another candidate with 
r3 ff . Obviously the present model is oversimplified to be 
applied to CeB 6 , especially because of the plane- wave 
conduction band with a single Fermi surface. We are now 



using more realistic band structures to derive intersite 
interactions. 

It is remarkable that the plane- wave model prefers r 5 „ 
octupole order if the intermediate states are domi- 
nant. The result is encouraging to understand the mi- 
croscopic origin of possible octupole order in phase IV 
of Ce x Lai_ x B 6 with x ~ 0.75. 14 " 16 ' It should be noted 
that the /° intermediate state does not favor r 5 „ at all. 
Therefore use of realistic energy band is again important 
to understand the microscopic mechanism. 

We now turn to CeB 2 C2 where the incommensurate 
magnetic order is observed. Since the Fermi surface is 
essentially an ellipsoid, our model may apply more di- 
rectly to this material. The lowest CEF states of CeB2C2 
is considered to be the m = ±1/2 doublet, which is split 
off from F8 by the tetragonal CEF. We have seen that 
the strong Kohn anomaly appears near the T point of 
the Brillouin zone. In view of the dHvA experiment of 
LaB 2 C 2 , 8 ) let us take 2k F = 0.86(27r/a g ) which gives 
5k = 0.14. The Kohn anomalies in our calculation arise 
at (<5"', 0,0) and (S"',S"',0) with S'" ~ 0.17 in the SCBZ. 
Mapping to the QCBZ, the relevant wavenumbers are 
(8 ,0,0), (S ,5 , 0), (S ,<5 ,<5 ) and their cyclic rota- 
tions. These values are in good correspondence with ex- 
perimental value of S = 0.16 already cited in §1. Actually 
the wavenumbers giving large intersite interaction form 
a ridge in the BZ, and the points quoted above are part 
of the ridge. We note that characteristic wavenumbers of 
RB2C2 compound, (1±<5, ±5, 0) in the QCBZ, 17 ) are close 
to intercepts of the sphere with radius 2/cf = 0.86(27r/a g ) 
at the BZ boundary. Thus we interpret incommensurate 
order found in RB2C2 in terms of the Kohn anomaly. 

In concluding this section we mention a theoretical 
study on the RKKY interaction between Pr nuclei, 18 ) 
where Kohn anomalies with a peak structure is found, 
although the authors do not remark on this feature. 

5. Summary 

We have calculated the RKKY interaction in the An- 
derson lattice with orbital degeneracy in the simple cubic 
lattice, and clarified the relation between the multipolar 
ordering structure and the Fermi surface. For intermedi- 
ate states, we have considered both 4/° and 4/ 2 Hund's- 
rule ground states as in the previous paper. Our simpli- 
fied model gives an incommensurate magnetic structure 
as the most stable structure, for which orbital angular 
momentum of 4/ electrons plays a crucial role. The prop- 
agation vector for incommensurate magnetic structure is 
determined by the Kohn anomaly at 2fcp which can be 
shifted by reciprocal lattice vectors. With one conduc- 
tion electron per cell, this anomaly can arise near the F 
point. If we take the Fermi wavenumber from the dHvA 
experiment of LaB 2 C2, the propagation vector is in good 
correspondence with that of incommensurate magnetic 
structure in CeB2C2- This result should be relevant in 
understanding the incommensurate magnetic structures 
observed in RB2C2 in general, but effects of anisotropy 
should be taken into account for more quantitative study. 

The anisotropy of CEF states influences the wavenum- 
ber dependence of interaction, i.e., the most favorable 
structure and the wavenumber of Kohn anomaly. The 
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propagation vector is unchanged whether the lowest CEF 
state is r 8 quartet or m = ±1/2 doublet for both inter- 
mediate states 4/° and 4/^. The most stable ordering 
structure in the case of 4/° intermediate state is the in- 
commensurate magnetism, while the 4/^ intermediate 
states give a staggered octupolar order. It is still an in- 
teresting open problem to clarify the relation between 
the multiplet structure of / shell other than Ce and the 
order parameters. 
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Appendix: Generalized Polarization Functions 

In this section, we present analytic results of gener- 
alized polarization functions. The radial parts of gener- 
alized polarization functions in the main text are given 

by 



/ 72*. (P*) / t2dt 
Jo z JO 



t 2 - 1' 



f 

Jo 



dz 

—Jh (pz)q(h,l 2 ;z). 



(A-l) 



where ji(x) is the spherical Bessel function of order I, and 
q(li,l 2 ; z)'s are tabulated in ref.9. Hereafter we consider 
even I3, since otherwise the integration about the solid 
angle il vanishes in eq.(8). 

From the analysis around z = 0, we find that 
z- 2 q(h,l 2 ;z) behaves as z ™™(h,h)+i ) an( j t ^ at in _ 
eludes terms like z - ' 3-1 sin(z) and z~' 3 cos(z). Therefore 
in reducing eq.(A-l) to Fourier sine and cosine trans- 
forms, the following three cases are considered: (i)Z 3 < 
min(/ 1 ,^ 2 ), (ii)?3 = l\ + l 2 , (iii)othcrwise. In the case (i), 
the calculation is straightforward once z~ 2 q(h, l 2 ; z) is 
divided into the sum of regular forms at z = 0, e.g. 



-^g(0, 0;z) = -^[cos(2z) 



7TZ^ 



2z l 



l]-^[sm(2z)~2z}. 



(A-2) 

In the case (ii), we transform the right hand side of 
eq.(A-l) as follows: 



f°° dz 
'Jo * 



ji 3 (z)q(h,h;z/p), 



(A-3) 



and consider the transform of ji s (z)/z 2 . This calcula- 
tion is almost the same as that for the case (i), because 



q(li,l 2 ;z)/z 2 includes terms like 



-l+i 



cos(z),z 'sin(z) 



with < I < (h + fo), and Si(z)/z 2 with the sine integral 
Si(z). In the case (hi), we combine the above two cases 
to obtain Xi\i 2 {p)- Relevant transforms are as follows: 



Jo 



dz 



sin (pz) 



z2 m+3 

(-l) m+1 
(2m + 2)! 



S in(z)-£(-iy 



y2 n+l 



n=0 



(2n+l)! 



(p+l) 2m+2 log|l + i| 
P 



-(P-1) 



2m+2 



log|l 



m+1 

E 

fc=0 



fc-1 



2ni f ..» \ p 



2fc-2i-l 



2k ) ^ 21 + 1 

l>0 



'2m + 2\ P 

k=0 x 7 i>0 



fc-1 



2fc-2i-l 



El LUl -f A \ 
\ 2k + 1 I 2^ ' 21 + 2 



(A-4) 



sin (pz) 



/ d *z 2 <<<+-> 



cos 



(z) -£(-!)' 



n=0 



(2n)! 



(_!)m+l n 



(2m +1)! [2 



(p+l) 2ro+1 lGg|l+-| 
P 



+ (p-l) 2m+1 log|l--| 

p 



k=0 



l>0 



21 + 1 



fc=0 



2m + 1 \ \ p 

21 + 2 
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l>0 



dz sin(pz) ( 4— I ^'( z ) 
\dz/ z 



/ dz sin(pz) I - 
Jo \ ( 

Jo 



dz cos(pz) 



d^ 211 - 1 Si(z) 
dz 



(A-6) 



dz cos(pz) [ 4— ) ^Kf) 
\dz J z 



(-1)" +1 
(2n + l) 2 



dz sin(pz) 



dz 



In 



Si(z) 



(A-7) 



where (™) is the binomial coefficient. 

In this way analytic result of x[\i 2 (p) is obtained, and 
some of the results are tabulated in Table A-l. In Ta- 
ble A-l we define functions f(x) and g(x) by 

1 + x 



f(x) = (x 2 - 1) log 
'°° dz 



9{x) 



1-x 
sin(xz)Si(z) 



(A-8) 



2 [L 2 {x) - L 2 {-x)] 0<|ar|<l 



[ 2± _ I 

= {i[L 2 (l/x)-L 2 (-l/x)} \x\>l 

where L 2 (x) is the Euler's dilogarithm, 19 ) 

log(l - z) 



L 2 (x) 



dz- 



(A-10) 



10 z 

We note that Xoo(p) 1S the Kohn function, and that 
Xu(p = 0) converges to the constant value, — 7r/4. 

Xh(p = 0) 



1 ^OO /-(X) 

= Hm / fc 2 dfc / fc' 2 dfc' 



Jo 



2j - /( £ fe') 



x dx 



k 2 - k' 



ji(kx)ji(k'x) 
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Table A-l. Polarization functions with orbital degeneracy 
2X; 3 ; (p)/n. The definitions of f(x) and g(x) arc given in the 
main text. This table gives results for max(£i, I2) < 2. 
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Here, we have used the following relation, 20 ) 
2a 2 Z" 00 

ji(aa;)ji(&a;)x 2 dx = 5 (a - b). 



2 j*l r 

71" Jo 



(A-ll) 



(A-12) 



This is related to the asymptotic behaviour of range func- 
tions in real space, 21 ) 



ii+i 



z- 4 q(h,l 2 ;z)~(-l) — 



2 z cos z — sin z 



(z — > 00). 
(A13) 

In general Xili 2 (p) is an oscillatory function for < p < 2, 
and are slowly varying for p > 2. We have checked the 
results in Fig. 3 or in Table A-l numerically by spherical 
Bcssel transforms. 22 ^ 
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